Phenotypic plasticity can facilitate adaptive evolution in gene regulatory circuits

BACKGROUND: Many important evolutionary adaptations originate in the modification of gene regulatory circuits to produce new gene activity phenotypes. How do evolving populations sift through an astronomical number of circuits to find circuits with new adaptive phenotypes? The answer may often involve phenotypic plasticity. Phenotypic plasticity allows a genotype to produce different - alternative - phenotypes after non-genetic perturbations that include gene expression noise, environmental change, or epigenetic modification. RESULTS: We here analyze a well-studied model of gene regulatory circuits. A circuit's genotype encodes the regulatory interactions among circuit genes, and its phenotype corresponds to a stable gene activity pattern the circuit forms. For this model, we study how genotypes are arranged in genotype space, where the distance between two genotypes reflects the number of regulatory mutations that set those genotypes apart. Specifically, we address whether this arrangement favors adaptive evolution mediated by plasticity. We find that plasticity facilitates the origin of genotypes that produce a new phenotype in response to non-genetic perturbations. We also find that selection can then stabilize the new phenotype genetically, allowing it to become a circuit's dominant gene expression phenotype. These are generic properties of the circuits we study here. CONCLUSIONS: Taken together, our observations suggest that phenotypic plasticity frequently facilitates the evolution of novel beneficial gene activity patterns in gene regulatory circuits

Protective effects of eicosapentaenoic acid plus hydroxytyrosol supplementation against white adipose tissue abnormalities in mice fed a high-fat diet

Objective: Obesity induced by high-fat diet (HFD) elicits white adipose tissue dysfunction. In this study, we have hypothesized that the metabolic modulator eicosapentaenoic acid (EPA) combined with the antioxidant hydroxytyrosol (HT) attenuates HFD-induced white adipose tissue (WAT) alterations. Methods: C57BL/6J mice were administered with a HFD (60% fat, 20% protein, 20% carbohydrates) or control diet (CD; 10% fat, 20% protein, 70% carbohydrates), with or without EPA (50 mg/kg/day), HT (5 mg/kg/day), or both for 12 weeks. Determinations in WAT include morphological parameters, EPA and docosahexaenoic acid content in phospholipids (gas chromatography), lipogenesis, oxidative stress (OS) and inflammation markers, and gene expression and activities of transcription factors, such as sterol regulatory element-binding protein-1c (SREBP-1c), peroxisome proliferator-activated receptor-gamma (PPAR-γ), nuclear factor kappa-light-chain-enhancer of activated B cells (NF-κB) (p65 subunit) and nuclear factor erythroid 2-related factor 2 (Nrf2) (quantitative polymerase chain reaction and enzyme linked immunosorbent assay). Results: HFD led to WAT hypertrophy in relation to PPAR-γ downregulation. WAT metabolic dysfunction was characterized by upregulation of lipogenic SREBP-1c system, mitochondrial energy metabolism depression, loss of the antioxidant Nrf2 signaling with OS enhancement, n-3 long-chain polyunsaturated fatty acids depletion and activation of the pro-inflammatory NF-κB system. EPA and HT co-supplementation diminished HFD-dependent effects additively, reaching values close or similar to controls. Conclusion: Data presented strengthen the importance of combined protocols such as EPA plus HT to attenuate metabolic-inflammatory states triggered by obesity

Boolean network model predicts cell cycle sequence of fission yeast

A Boolean network model of the cell-cycle regulatory network of fission yeast (Schizosaccharomyces Pombe) is constructed solely on the basis of the known biochemical interaction topology. Simulating the model in the computer, faithfully reproduces the known sequence of regulatory activity patterns along the cell cycle of the living cell. Contrary to existing differential equation models, no parameters enter the model except the structure of the regulatory circuitry. The dynamical properties of the model indicate that the biological dynamical sequence is robustly implemented in the regulatory network, with the biological stationary state G1 corresponding to the dominant attractor in state space, and with the biological regulatory sequence being a strongly attractive trajectory. Comparing the fission yeast cell-cycle model to a similar model of the corresponding network in S. cerevisiae, a remarkable difference in circuitry, as well as dynamics is observed. While the latter operates in a strongly damped mode, driven by external excitation, the S. pombe network represents an auto-excited system with external damping.Comment: 10 pages, 3 figure

Network model of immune responses reveals key effectors to single and co-infection dynamics by a respiratory bacterium and a gastrointestinal helminth

Co-infections alter the host immune response but how the systemic and local processes at the site of infection interact is still unclear. The majority of studies on co-infections concentrate on one of the infecting species, an immune function or group of cells and often focus on the initial phase of the infection. Here, we used a combination of experiments and mathematical modelling to investigate the network of immune responses against single and co-infections with the respiratory bacterium Bordetella bronchiseptica and the gastrointestinal helminth Trichostrongylus retortaeformis. Our goal was to identify representative mediators and functions that could capture the essence of the host immune response as a whole, and to assess how their relative contribution dynamically changed over time and between single and co-infected individuals. Network-based discrete dynamic models of single infections were built using current knowledge of bacterial and helminth immunology; the two single infection models were combined into a co-infection model that was then verified by our empirical findings. Simulations showed that a T helper cell mediated antibody and neutrophil response led to phagocytosis and clearance of B. bronchiseptica from the lungs. This was consistent in single and co-infection with no significant delay induced by the helminth. In contrast, T. retortaeformis intensity decreased faster when co-infected with the bacterium. Simulations suggested that the robust recruitment of neutrophils in the co-infection, added to the activation of IgG and eosinophil driven reduction of larvae, which also played an important role in single infection, contributed to this fast clearance. Perturbation analysis of the models, through the knockout of individual nodes (immune cells), identified the cells critical to parasite persistence and clearance both in single and co-infections. Our integrated approach captured the within-host immuno-dynamics of bacteria-helminth infection and identified key components that can be crucial for explaining individual variability between single and co-infections in natural populations

Most Networks in Wagner's Model Are Cycling

In this paper we study a model of gene networks introduced by Andreas Wagner in the 1990s that has been used extensively to study the evolution of mutational robustness. We investigate a range of model features and parameters and evaluate the extent to which they influence the probability that a random gene network will produce a fixed point steady state expression pattern. There are many different types of models used in the literature, (discrete/continuous, sparse/dense, small/large network) and we attempt to put some order into this diversity, motivated by the fact that many properties are qualitatively the same in all the models. Our main result is that random networks in all models give rise to cyclic behavior more often than fixed points. And although periodic orbits seem to dominate network dynamics, they are usually considered unstable and not allowed to survive in previous evolutionary studies. Defining stability as the probability of fixed points, we show that the stability distribution of these networks is highly robust to changes in its parameters. We also find sparser networks to be more stable, which may help to explain why they seem to be favored by evolution. We have unified several disconnected previous studies of this class of models under the framework of stability, in a way that had not been systematically explored before

The logic of the floral transition: reverse-engineering the switch controlling the identity of lateral organs

Much laboratory work has been carried out to determine the gene regulatory network (GRN) that results in plant cells becoming flowers instead of leaves. However, this also involves the spatial distribution of different cell types, and poses the question of whether alternative networks could produce the same set of observed results. This issue has been addressed here through a survey of the published intercellular distribution of expressed regulatory genes and techniques both developed and applied to Boolean network models. This has uncovered a large number of models which are compatible with the currently available data. An exhaustive exploration had some success but proved to be unfeasible due to the massive number of alternative models, so genetic programming algorithms have also been employed. This approach allows exploration on the basis of both data-fitting criteria and parsimony of the regulatory processes, ruling out biologically unrealistic mechanisms. One of the conclusions is that, despite the multiplicity of acceptable models, an overall structure dominates, with differences mostly in alternative fine-grained regulatory interactions. The overall structure confirms the known interactions, including some that were not present in the training set, showing that current data are sufficient to determine the overall structure of the GRN. The model stresses the importance of relative spatial location, through explicit references to this aspect. This approach also provides a quantitative indication of how likely some regulatory interactions might be, and can be applied to the study of other developmental transitions

Flower Development as an Interplay between Dynamical Physical Fields and Genetic Networks

In this paper we propose a model to describe the mechanisms by which undifferentiated cells attain gene configurations underlying cell fate determination during morphogenesis. Despite the complicated mechanisms that surely intervene in this process, it is clear that the fundamental fact is that cells obtain spatial and temporal information that bias their destiny. Our main hypothesis assumes that there is at least one macroscopic field that breaks the symmetry of space at a given time. This field provides the information required for the process of cell differentiation to occur by being dynamically coupled to a signal transduction mechanism that, in turn, acts directly upon the gene regulatory network (GRN) underlying cell-fate decisions within cells. We illustrate and test our proposal with a GRN model grounded on experimental data for cell fate specification during organ formation in early Arabidopsis thaliana flower development. We show that our model is able to recover the multigene configurations characteristic of sepal, petal, stamen and carpel primordial cells arranged in concentric rings, in a similar pattern to that observed during actual floral organ determination. Such pattern is robust to alterations of the model parameters and simulated failures predict altered spatio-temporal patterns that mimic those described for several mutants. Furthermore, simulated alterations in the physical fields predict a pattern equivalent to that found in Lacandonia schismatica, the only flowering species with central stamens surrounded by carpels

Can computational efficiency alone drive the evolution of modularity in neural networks?

Some biologists have abandoned the idea that computational efficiency in processing multipart tasks or input sets alone drives the evolution of modularity in biological networks. A recent study confirmed that small modular (neural) networks are relatively computationally-inefficient but large modular networks are slightly more efficient than non-modular ones. The present study determines whether these efficiency advantages with network size can drive the evolution of modularity in networks whose connective architecture can evolve. The answer is no, but the reason why is interesting. All simulations (run in a wide variety of parameter states) involving gradualistic connective evolution end in non-modular local attractors. Thus while a high performance modular attractor exists, such regions cannot be reached by gradualistic evolution. Non-gradualistic evolutionary simulations in which multi-modularity is obtained through duplication of existing architecture appear viable. Fundamentally, this study indicates that computational efficiency alone does not drive the evolution of modularity, even in large biological networks, but it may still be a viable mechanism when networks evolve by non-gradualistic means